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Credible control and overall stabilization of closed-loop nonlinear systems 
presented in a port-controlled hamiltonian structure (PCH-D) were forever 
the development of the simultaneous interconnection and damping 
assignment passivity-based control (SIDA-PBC) method. The robustness 
and reliability of the method called into question against noise and certainly 
modelling errors. Indeed, a new scheme has been presented to control a 
doubly fed induction generator (DFIG) based on the energy form and 


Keywords: exploiting the electrical parameters of the closed loop system. The results 
DFIG obtained provide a new technique and implement more freedom when 

a designing the diagram of the advanced controller during the production of 
Passivity 


active power. The contribution of this paper is researching on the advanced 
PCH model nonlinear methodes used, many simulations were carried out in simulation 
SIDA-PBC control using the MATLAB/Simulink environment under important operating 
conditions, allowing to demonstrate the feasibility of the proposed method 
and verify the performance considering the robustness. 
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1. INTRODUCTION 

Passive-based control (PBC) has become a research tool for nonlinear control, mainly because of its 
ease of application to physical systems (mechanical, electrical, and electromechanical) [1], [2]. In recent 
years, assignment-passivity-based interconnection and control (SIDA-PBC) has been developed as a flexible 
and general nonlinear system controller design method, and introduced tools that affect interconnection and 
internal energy loss considerations [2]. This approach enables a robust controller with a clear physical 
interpretation of the system's connection to its environment [3]-[6]. In particular, the total energy of the 
closed loop is the difference between the system energy and the energy provided by the controller [7]-[10]. 
Furthermore, power generation is an interesting research area for nonlinear control. They are especially used 
by doubly fed induction generator (DFIG). 

This research analyzes the use of SIDA-PBC as a method to control the power provided by this 
entrainement. The result is a new approach to designing controllers for these workouts, focusing on energy 
characteristics such as balance and fitness to achieve desired goals. Therefore, a power controller was 
obtained from the new scheme and simulated using the DFIG model in the MATLAB Simulink environment. 
To apply the SIDA-PBC method to generator induction, a connection control Hamiltonian (PCH) model is 
derived for the entire electromechanical system [2], [11]. The SIDA-PBC controller decouples the input of 
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the induction generator, exploiting the internal characteristics of the system to achieve the desired goal 


without intermediate steps [12]—[15]. 


2. DFIG: PCH MODELISATION WITH SIDA-PBC 
2.1. Machine modelisation 
Willingly, the PARK coordinate system linked to the rotating field of the DFIG is written as (1) [16]. 
Vas =R slas gee esas 
Vas = Rslqs pej wsPas 
ia — (Ws — ©) P gr (1) 
Vor = Rylqr + ae + (ws — w) Par 


d 
JZ = M(larlqs — laslqr) — Cr — Cfo 


Var = Rrlar + 


and flux equations can be written as (2). 


as = Vas—Rslas + wsLslqs + wsMlgr 
Pqs = Vgs —Rslqs — WsLslas — wsMlar 


; 2 
Par = Var—Rylar + (ws — @)Lylgr + (Ws — w)Mlgs 2) 
Do, > = Vor Ry gr — (ws ai w)Lrlar = (ws a w)Mlas 
where: 
, = VŅ-R;hl; — wsLsJzl; — wsMjJzl, (3) 
ò, rR lal, — (ws ~ @)MJols = (ws 7 w)L Joly (4) 
and, 
d 
Jorie op = LsrlsJaly — C, — Cro (5) 
p P I I = 
with: ®, = l ©, = -fochs = L= ehes = al pe [5 abd: = h a] 


The state variables are: 


x= [o ,®, T Jorreo| = eee Xm IF 


with: x7 = [®,”, o,."| , Xm = Jprigw: are the electrical and Mechanical state variables respectively. The 
passive energy equation as in (6). 


H(x) = 5x01 yo + x2, (6) 
2JDFIG 
wp. 7 _ [Esk E 
With: L = Mb Lh 


For the previous equation, the derivatives of electrical and mechanical state variables as in (7). 


ð 
se 0. (7) 
=f H 1 
— =J x — 
OXm DFIG“m 


Lastly, the matrices control written as: 


—wsLsJ2 —w;M]J2 02x1 Rsl2 02x2 02x1 
I(x) = |—osLsJz —(ws — w)LsJz MJzls|, R(x) = |02x2 Rrl2 02x1 
01x2 MII J, 0 Oix2 Oix2 G 
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h 02x2 02x1 


g(x) = [os h 02x1 


| u= [v7 ye T1” 
01x2 01x2 1 


wwie ohh = [0 OOA RORE 


Conclusively, the mathematical control is: 


—wsLsJ2 —w;MJ2 02x1 Rsl2 02x2 02x1 
X= ||-osMJ2 —(ws— w)LsJz MJzls|— |02x2 Retz 02x1ı||VH 
Oixz MII), 0 Oix2 Oi1x2 Ty 
h 02x2 02x1] [Kr 

+102%2 hk — Ozxı| |V 
0ix2 01x2 1 T, 
h 02x2 02x1 
y= 02x2 h 02x1 VH 
Oix2 01x2 1 


2.2. SIDA-PBC methode 
The energy equation as in (8) [3], [12], [17]-[22]. 


Ha(x) = 527 Pak, Py = PZ > 0 (8) 
SIDA-PBC methode is solving Fy(x) and countering Fy(x) + FẸ < 0. 
V; 02x2 
VU, w) — R]ƏH + ]02x1) +| Lp |V = F(x) Paf (9) 
Cem 01x2 
To simplify the solution, we offer P4 diagonal. 
PS 02x2 02x1 Falx) Fz)  Ozx1 
Pa = |02x2 Ps Ozx1| > 0, Fa(x) = |F) Fz2(x) Fz3(x) 
Oix2 Oix2 Po FLC) FL) F33(x) 
Where: F(x) is suit P4. 
By the (9) and the equilibria (7), we have: 
—=(wsLsJz + Relay — wsLsJ2l, = (LFP; + MF,2P,)Î; + (MFR + Ly Fi2P, Jl, (10) 
For everything /,, [., there is a unique solution given by: 
1 Lr 
Fii = -4 (osz +ZR,h) (11) 
(12) 


M 
Fiz = putsh 


With: u = L,L, — M? > 0 


The calcul of F32 of Ja(x) is possible set to zero. 
The choice of F3,(x) be easy by having this precedent term equal to zero. Of course, it appears in 


Fy(x) + FẸ (x), and F32(x) is not critical because his contribution to the F(x) + Fj (x) perhaps centered by 
an appropriate choice of F,3(x) that given the presence of the command [22]-[24]. 


[FP F3 (x)P-]LI + F33(x)PaJprig® = MI Jz, — Cro = M- hli- Cr@ (13) 
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From the term @ we obtain: 
Fs; = ———— (14) 


PwIDFIG 


While the remaining equations can be arranged in ways: 


PF34(x)] _ u| eE ) 
=0 15 
(1 BERO "FRL (15) 
7 — MEJ Mh 
With: L = |M al 


By adding to (10) a vector G(x) E€ R* 


ar (, [BF (2) BE _ 
w (x aed =M al Š Ga) =A ue 
Or: i. G(x) = 0 


and fixing the front of (16)=0 we have: 


PSF 33(x) M L;J2lų + MJols s |+ +e c(X) (17) 
P,Fz2(x) u \—LrJz2ls — MJol, Gp (x) 
and by introducing the therm: 
Ge(X)] _ 7-1 
Go| = GC) (18) 


With: 
G_(x), Gp(x) E R? 


As mentioned above, in order to satisfy the antisymmetry condition, it is necessary to find a solution 
with F3,(x) = cst, which can be easily achieved by choosing [16], [23]-[25]. 


M? > 
G(x) = Ses (19) 
With this equation, we result: 


ay I, + MGp 
G(x)=| # (20) 
=gh! + L,Gp 


and, in order to ensure if G(x) = 0 , we fix: 
M2. > 
Gp(x) = — Sale (21) 


Lastely, by replacing (19) and (21) in (17) we have: 


M * * M * 
F3, = z2 Ms + L, Ir) = pee (22) 
M M 
F32 = — 5 plallsls + MI,) = me les (23) 


In the next step, the item F4(x), is selected, which is directly affected by the control measure [16], [24]-[26]. 
To satisfy the symmetry constraints of Jq, and simplify the state, we choose: 


kr 
Foy = —Fy2 , Fo3(x) = —F32(x) , F22 Fp <0. 
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With: 
2LrRs M š 
= p R 02x2 pg )2®r 
kr 
Fa(x) + FF (œ) =| 02x2 = z2 02x1 (24) 
V” ox ef 
Psh Da 01x2 PwIJDFIG 


Analysis to determine this is easy to complete Fy(x) + Fj (x) < 0 if and only if the parameters P, and P, 
satisfied: 


JorigM? 
P, > {2S p4] |p 
s (a G 


The equations of tensions of the DFIG in closed loop is written as: 
V = Rl, + (ws ant w)J2(MI1, + L,I,) + K,(LsI, + MI,) > K, (LI, + MI,) + KoJ2 Ps6 (25) 
With: 


PoloricM M? 
K, > 0,K, =————_ > 0 et K, > ——_——_|9; |? K, 
r w Pu et Ks 4GL,L,n' rl w 


The global system control diagram is shown in Figure 1: 


SIDA-PBC Control 
of Power 


Figure 1. SIDA-PBC control scheme 


3. DISCUSS THE RESULTS 

Using the MATLAB environment simulation and with the parameters of the doubly fed induction 
generator (DFIG) mentioned in Table 1 in the appendix; we therefore obtained in Figure 2 the tests results of 
the track and regulation for the SIDA-PBC control. Figure 3 presents the results of the robustness tests, 
Figure 3(a) results of 30% variation of inductance and Figure 3(b) results of 30% variation of resistance. 


Table 1. DFIG parameters [2] 


Parameter Value 
nominal power 100 Kw 
Stator resistance 0.40 
Rotor resistance 0.10 
Stator inductance 0.07 H 
Rotor inductance 0.02 H 
Mutual inductance 0.03 H 
Moment of inercia 0.5 Kg.m? 


coefficient of friction 0.002 N.m.s/rad 
number of pole pairs 2 
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Figure 2. Tests results of the track and regulation 
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Figure 3. Robustness tests (a) results of 30% variation in inductance and (b) results of 30% change in 
resistance 


4. CONSULTATION 

Evaluation of power controller performance by nonlinear control of the generator performance by 
the controller based on switching and damping (SIDA-PBC) considerations. We ran a series of 
MATLAB/Simulink simulations. In a tracking test, the controlled variable successfully follows the target 
path without exceeding the active and reactive powers and without steady-state static errors. The 3-phase 
stator current draw is low at the start and sinusoidal at steady state. The rotor current frequency is 17.5 Hz 
thanks to a good choice of PWM control strategy in the inverter. Proportional electromagnetic torque and 
active power have the same appearance at constant speed. 

In simulations tests, we notice very good sensitivity to disturbances caused by abrupt setpoint 
changes from -20 KW to -80 KW, active and reactive powers are unaffected. The three-phase rotor current is 
sinusoidal. Electromagnetic torque and active power output are not affected. Robustness tests show that a 
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+50% variation in resistance (Rs, Rr) and +30% in inductance (Ls, Lr) has very little effect on response time 
and osciations amplitudes. 


5. CONCLUSION 

The presented research focuses on active and reactive powers advanced control of a doubly fed 
induction generator by a new structure of passivity. Robustness results show a very good sensitivity. Active 
and reactive powers always stay within tolerance for disturbances this means that the variation of the 
parameters (resistances and inductances rotorique and statorique) little impact on response time. It has been 
shown that very good results can be obtained thanks to the UPS sine-triangle PWM control strategy. 
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